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Abstract 

The recent development of a micromagnetic simulation methodology — suitable for multiphase magnetic 
nanocomposites — permits the computation of the magnetic microstructure and of the associated magnetic small-angle 
neutron scattering (SANS) cross section of these materials. In this review article we summarize results on the micro- 
magnetic simulation of magnetic SANS from two-phase nanocomposites. The decisive advantage of this approach 
resides in the possibility to srutinize the individual magnetization Fourier contributions to the total magnetic SANS 
cross section, rather than their sum, which is generally obtained from experiment. The procedure furnishes unique 
and fundamental information regarding magnetic neutron scattering from nanomagnets. 
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1. Introduction 

Small-angle neutron scattering (SANS) is one of the 
most important techniques for microstructure determi- 
nation in soft and hard condensed matter, materials sci- 
ence, and in physical chemistry. Since most SANS stud- 
ies focus on nuclear rather than magnetic scattering, it 
is not surprising that the theoretical concepts behind 
nuclear SANS are rather well developped ljj-|9|]. By 
contrast, the understanding of magnetic SANS is still 
at its beginning, although magnetic SANS has previ- 
ously demonstrated great potential for resolving the spin 
structures of various magnetic materials. For instance 
(in the last decade), magnetic SANS has been employed 
for stud ying th e microstructures of magnetic nanocom- 
posites lll0l - ll9ll . amorphous alloys 120 42311 . and of ele- 
mental nanocrystalline bulk ferromagnet s 124 - 13511 . the 
process of dynamic nuclear polarization 11361] . imaging 
of the flux-line lattice in superconductors 0371 13811 . pre- 
cipitates in steels 13911 . nanocrystalline rare-earth metals 
with random paramagnetic susceptibility 14011 . fractal 
magnetic domain structures in NdFeB permanent mag- 
nets 114 111 , spin structures of ferro fluids, nanoparticles, 
and nanowires 142145211 . magnetostriction in FeGa alloys 
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115 311 . electric-field-induced magnetization in multifer- 
roics f54ll. magnetization reversal in magnetic recording 
media [55] and exchange-bias materials [56], and chiral 
and skyrmion-like structures in single crystals ll57l - l59ll . 

Necessary prerequisite for the quantitative analysis 
of elastic magnetic SANS data is the knowledge of 
the Fourier components of the static magnetization vec- 
tor field M(r) of the sample under study. The theory 
of micromagnetics [60-62] provides the proper frame- 
work for the computation of M(r). However, the solu- 
tion of Brown's equations of micromagnetics amounts 
to the solution of a set of nonlinear partial differen- 
tial equations with complex boundary conditions, a task 
which cannot be done analytically for most practically 
relevant problems. Therefore, closed-form expressions 
for the ensuing so-called spin-misalignment scattering 
cross section are limited to the approach-to-saturation 
regime 113 ll 16311 . in which the micromagnetic equations 
can be linearized. 

In this review article we summarize our recent work, 
in which we have used numerical micromagnetics for 
the computation of the magnetic SANS cross section of 
two-phase magnetic nanocomposites. The use of nu- 
merical techniques allows us to solve the underlying 
equations rigorously, without resorting to the high-field 
approximation (saturation regime). This approach pro- 
vides insights into the fundamentals of magnetic SANS. 
The micromagnetic simulations are adapted to the mi- 
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crostructure of a two-phase nanocomposite from the 
NANOPERM family of alloys Q. 

The paper is organized as follows: in Sec. 2 we 
provide the details of our micromagnetic methodology, 
Sec. 3 discusses the magnetic SANS cross sections for 
the two most commonly used scattering geometries, 
Sec. 4 presents and discusses the results of the micro- 
magnetic simulations for the magnetic SANS cross sec- 
tion, and in Sec. 5 we summarize the main findings of 
this study. 

2. Micromagnetic background 

Micromagnetism is a mesoscopic phenomenological 
theory designed to compute the equilibrium magnetiza- 
tion state of an arbitrarily shaped ferromagnetic body, 
when the applied field, the geometry of the ferromagnet 
and all materials parameters are known 1 60-62]. In or- 
der to find the equilibrium magnetization configuration 
Meq(r), the total magnetic free energy of a ferromag- 
net should be considered as a functional of its magne- 
tization state, E tot = Z? tot [M(r)]. The state which de- 
livers a (local) minimum to this functional corresponds 
to the required equilibrium magnetization configuration, 
so that the problem amounts to the minimization of the 
total energy functional. In the most common case E tol 
contains contributions from the energy due to an exter- 
nal field, exchange, anisotropy and magnetodipolar in- 
teraction energies. Due to the nonlocal nature of the 
magnetodipolar interaction, almost all practically inter- 
esting problems can not be treated analytically, so that 
numerical minimization of £' tot [M(r)] should be carried 
out. In the contemporary research landscape, numerical 
micromagnetics is a large and still continuously expand- 
ing field. Recent reviews on the micromagnetic state of 
the art can be found in the handbook Ref . B65I1 . In this 
article we briefly discuss only those methodical aspects 
of numerical micromagnetics which are important for 
simulations of nanocomposite materials. 

First, we would like to emphasize that such materials 
are one of the most complicated objects from the point 
of view of numerical simulations. The main difficulty is 
that they consist of at least two phases, and the bound- 
aries between these phases are complicated curved sur- 
faces; a typical example is a hard-soft nanocompos- 
ite consisting of magnetically hard (i.e., having a large 
magnetocystalline anisotropy) crystal grains surrounded 
by a magnetically soft matrix. Such a system is very 
difficult to simulate for the following reasons. The 
majority of modern numerical micromagnetic methods 
can be subdivided into two classes, the so-called finite- 
difference and finite-element methods (FDM and FEM) 



Il65ll . In FDM the system under study is discretized into 
a regular translationally invariant (usually rectangular) 
grid. Such a discretization allows, first, the evaluation of 
the exchange field by simple finite-difference formulas, 
which are the finite-difference approximations for the 
corresponding second-order differential operator acting 
on the magnetization field M(r) (see below). Second, 
the translational invariance of a FDM grid enables the 
usage of the fast Fourier transformation (FFT) technique 
for the computation of the long-range magnetodipolar 
interaction field and energy. For a system discretized 
into N cells, the FFT technique reduces the operation 
count for this energy from ~ A^ 2 (for a direct summa- 
tion) to ~ N log N. However, a serious disadvantage 
of a regular grid is a pure approximation for arbitrarily 
curved surfaces and boundaries. This is an important 
drawback for simulations of magnetic nanocomposites, 
because the adequate representation of the interphase 
boundaries for the accurate evaluation of associated ex- 
change and magnetodipolar interactions between differ- 
ent phases is crucially important. 

The second group of numerical methods widely used 
in micromagnetics — finite-element methods (FEM) — 
employ the discretization of the system under study 
into arbitrarily shaped tetrahedrons. The flexibility of 
this discretization type allows one to represent curved 
boundaries (including those between magnetically hard 
inclusions and the soft magnetic matrix) with any de- 
sired accuracy. However, the price to pay for this flex- 
ibility is high. First, computation of the exchange field 
requires now complicated methods designed for the ac- 
curate representation of second-order differential opera- 
tors on irregular lattices. Second (and most important), 
it is no longer possible to use FFT for the magnetodipo- 
lar field evaluation. For this reason, highly sophisticated 
methods for the computation of this field are used in 
FEM simulations. These methods, which are based on 
the decomposition of magnetic potentials inside the fer- 
romagnet and in the outer space, and the subsequent so- 
lution of the corresponding Poisson equations for these 
potentials on irregular grids [65] require a high pro- 
gramming effort and result in iterative algorithms for 
the evaluation of the dipolar field for a given magneti- 
zation configuration (in contrast to the FFT technique). 

Another important limitation of finite-element meth- 
ods is that they can only be employed in simulations 
with open boundary conditions (OBC), so that periodic 
boundary conditions (or PBC, routinely applied in simu- 
lations of extended thin films and bulk materials in order 
to eliminate strong finite-size effects) can not be used. 
The impossibility to apply PBC is a serious disadvan- 
tage in simulations of SANS experiments on nanocom- 
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posites, whereby the scattering intensity is sensitive to 
magnetization fluctuations in the bulk. Artificial sur- 
face demagnetizing effects arising in simulations with 
OBC might be very significant in this case, due to a 
relatively small simulation volume affordable even for 
modern computers. In addition, the suppression of these 
effects is especially important for nanocomposites con- 
taining a soft magnetic phase. 

Another undesirable feature of a tetrahedron mesh is 
that hard magnetic grains must also be discretized into 
tetrahedrons, although in many cases the magnetization 
within a single grain is nearly homogeneous. This leads 
to a significant increase of the total number of finite ele- 
ments required, resulting in a corresponding increase of 
the computation time; we refer the reader to Ref. ll66ll 
for the discussion of this problem. 

Due to all the reasons explained above, numerical 
micromagnetic simulations of SANS experiments on 
nanocomposites are very rare [67, 68]. Corresponding 
full-scale simulations of SANS measurements on a two- 
phase system have been reported, up to our knowledge, 
only in Ref. [67], where the magnetization configuration 
of a longitudinal magnetic recording media film was 
modeled. Based on the experimental characterization 
of this material, the authors of l67ll have built a two- 
phase model for this system, where each magnetic grain 
consisted of a hard magnetic grain core and an essen- 
tially paramagnetic grain shell, having a very high sus- 
ceptibility. The OOMMF code employing the standard 
FDM has been used [69], so that a very fine discretiza- 
tion (0.3 x 0.3 x 0.3 nm 3 cells) had to be applied in order 
to reproduce the spherical shape of grain cores with a 
required accuracy. For this reason, only a rather limited 
number of grains (~ 50) could be simulated. In addi- 
tion, the exchange interaction both between the grains 
and within the soft magnetic matrix (represented by the 
merging grain shells) was neglected. Still, using sev- 
eral adjustable parameters, a satisfactory agreement of 
the simulated SANS intensity profile with experimental 
data was achieved. 

The brief overview of the methodical problems pre- 
sented here clearly shows that both a qualitative im- 
provement of the micromagnetic simulation methodol- 
ogy and extensive numerical studies devoted to SANS 
experiments are highly desirable. 

2.1. New micromagnetic algorithm: mesh generation 
and its regular representation 
For the reasons explained in the previous section and 
in order to perform accurate and efficient simulations of 
two-phase nanocomposites, we need to generate a poly- 
hedron mesh with the following properties: (a) it should 



allow to represent each hard magnetic crystallite as a 
single finite element (because the magnetization inside 
such a crystallite is essentially homogeneous), (b) the 
mesh should allow for an arbitrarily fine discretization 
of the soft magnetic matrix in-between the hard grains 
(to account for the possible large variations of the mag- 
netization direction between the hard grains), and (c) the 
shape of the meshing polyhedrons should be as close as 
possible to spherical, in order to ensure a good quality 
of a spherical dipolar approximation for the calculation 
of the magnetodipolar interaction energy, even for the 
nearest neighboring mesh elements. 

A mesh consisting of polyhedrons satisfying all these 
requirements can be generated using two kinds of meth- 
ods. First, there exist various modifications of a purely 
geometrical algorithm designed to obtain a random 
close packing of hard spheres [70]. In these algorithms, 
the initial distribution of sphere centers is completely 
random. Then, at each step the worst overlap between 
two spheres is eliminated by pushing these spheres apart 
along the line connecting their centers. This procedure 
usually introduces new overlaps; however, these over- 
laps are usually smaller and are eliminated during the 
next steps, so that on the average the packing quality 
improves (the largest overlap present in the system de- 
creases). The algorithm is robust and produces a ran- 
dom close packing of nonoverlapping spheres with any 
desired accuracy (see Ref. tp70h for further details). Un- 
fortunately, the computation time for this method in- 
creases with the number of elements N as ~ N 2 , so that 
the maximal number of spheres which can be positioned 
within a reasonable computation time is N ~ 10 4 . 

Therefore, in order to generate a mesh with a much 
larger number of finite elements (N > 10 5 ), we have de- 
veloped a "physical" method, where we model a system 
of spheres interacting via a short-range repulsive poten- 
tial: 



Zf dj; - (fi + rf) 
A pot exp ' 



n. 



ec 



(l) 



Here, the constant A pot determines the value of our po- 
tential when the distance d/j between the centers of in- 
teracting spheres is equal to the sum of their radii r, and 
rj (to ensure small overlaps in the final configuration, 
it should be A pot » 1; in a typical case A pot = 10). 
The parameter rd ec defines the decay radius of the po- 
tential. Again, at the beginning of iterations, sphere 
centers are positioned randomly. Then, we move the 
spheres according to the purely dissipative (i.e., neglect- 
ing the inertial term) equation of motion resulting from 
the forces obtained from the potential Eq. (Q}. The time 
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step for the integration of this equation is adjusted to 
ensure decrease of the total system energy after each 
step. Due to the repulsive nature of the potential Eq. (Q]), 
this procedure leads also to the decrease of overlaps of 
the spheres. To achieve the desired result, we move the 
spheres until their maximum overlap does not exceed 
some prescribed small value (we have found that for 
our purposes the remaining overlap (r, + rj)/dij > 0.95 
is good enough). The algorithm may be refined further 
to increase its efficiency; in particular, one might de- 
crease the decay radius of the potential r^ ec , thus mak- 
ing the potential "harder", when the overlapping be- 
tween spheres decreases during the sphere motion. We 
also note that due to the random spatial arrangement of 
spheres obtained in this way, we avoid possible artifacts 
caused by the regular placement of finite elements. 

After the spheres have been positioned using one of 
the two algorithms described above, their centers are 
used as location points of magnetic dipoles. To compute 
the magnitudes of these dipoles, we have to determine 
the volume of each corresponding mesh element, which 
is in fact a polyhedron. This determination is made via a 
regular grid representation procedure that should satisfy 
the following requirements. First, we should conserve 
the total sample volume. Second, the interface between 
neighboring mesh elements should be fiat as far as pos- 
sible (apart from geometrical reasons). The last require- 
ment is also supported by electron microscopy images 
of various polycrystalline magnets (e.g., J7I, 72]). 

In order to satisfy both these requirements, we used 
the following method: the sample is divided into cubi- 
cal cells which side is much smaller (usually about four 
times smaller) than the size of a finite element (polyhe- 
dron) of our disordered mesh used to discretize the soft 
phase. For every cubical cell (j, k, I), we calculate the 
distances As'. , , between the center of this cell and the 
centers of neighboring polyhedrons (labeled by i). The 
function 



min (; , [(As' y) 2 -^] 



(2) 



indicates to which polyhedron (with radius /?,-) we at- 
tribute the current (J, k, I) cube. The sum of cube vol- 
umes ascribed to the given polyhedron is taken as its 
volume. As a result of this procedure, the distribu- 
tion of mesh-element volumes for both magnetic phases 
demonstrates a nearly Gaussian behavior. To obtain the 
magnitude of the dipolar moment assigned to each poly- 
hedron, its volume is multiplied by the saturation mag- 
netization of the material inside which the polyhedron is 
located (we remind that nanocomposites consist of ma- 
terials with different magnetizations). 




Figure 1 : Schematical representation of the mesh-generation method: 
spheres in the left image indicate the distribution of magnetic 
dipoles (blue — hard magnetic phase, yellow-orange-red — soft mag- 
netic phase). The corresponding regular grid representation on the 
right is used for the mesh-element volume determination (see Sec. 12.11 
for details). Note that the actual "sample", which is used for the mi- 
cromagnetic simulations, is featured in Fig. [3] 



This method also allows for a very efficient calcu- 
lation of the Fourier components of the magnetization 
(see Eqs. ( TTTb and ( TT8b below) for a disordered system, 
using FFT on the already composed regular grid. 

Summarizing, the whole algorithm can be viewed as 
a method to discretize a sample into polyhedrons hav- 
ing nearly spherical shape (see Fig. Q]). This is due to 
the fact that polyhedrons "inherit" the spatial structure 
obtained by the positioning of closely packed spheres. 
The fact that the shape of the volume which is occupied 
by each magnetic moment is nearly spherical allows us 
to use the spherical dipolar approximation (equivalent 
to the point dipole approximation) for the evaluation of 
the magnetodipolar interaction between the moments. 

Finally, we point out that both algorithms allow for 
the usage of polyhedrons with different sizes, if we need 
different meshing on different system locations. 

2.2. New micromagnetic algorithm: energy contribu- 
tions 

In our micromagnetic simulations we take into ac- 
count all four standard contributions to the total mag- 
netic free energy listed above: energy in the exter- 
nal magnetic field, energy of the magnetocrystalline 
anisotropy, exchange stiffness and magnetodipolar in- 
teraction energies. 

2.2.1. External field and magnetocrystalline anisotropy 
energies 

The system energy due to the presence of an external 
magnetic field and the energy of the magnetocrystalline 
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anisotropy (which can be uniaxial and/or cubic) are cal- 
culated in our model in the standard way, namely 



£ext = -^//,H, 
i=i 

N 

£™ = -^^ n V,(m,n,) 2 , 



(3) 



(4) 



where H is the external field, //, = //(r,) and V, are the 
magnetic moment and the volume of the z'th finite ele- 
ment (polyhedron), and m, denotes the unit magnetiza- 
tion vector. Both the anisotropy constants Kj and the di- 
rections of the anisotropy axes n, can be site-dependent, 
as required for a polycrystalline nanocomposite mate- 
rial. For the cubic anisotropy case, the symbols m^ x < etc. 
represent the components of unit magnetization vectors 
in the local coordinate system that is attached to the cu- 
bic anisotropy axes. 

2.2.2. Exchange energy 

The evaluation of this energy contribution in our 
model requires a much more sophisiticated approach 
than in the standard FDM, because the continuous in- 
tegral version of this energy contains the magnetization 
gradients, 

£exch = J A(r) [(Vm,) 2 + (Vm v ) 2 + (Vm e ) 2 ] dV,(6) 

where A denotes the exchange-stiffness parameter and 
V is the sample volume. Finding an approximation to 
Eq. (0 for a disordered system is a highly nontrivial 
task. 

We remind that for a regular cubic grid with a cell size 
a (and cell volume AV = a 3 ), it can be shown rigorously 
(see the detailed proof in Ref. [73]) that the integral in 
Eq. (0 can be approximated as the sum 



^exch 



jZZ 

i=l jcn.n.(i) 



2A U AV i N 
—L-( m mj). (7) 



Here, Ay denotes the exchange-stiffness constant be- 
tween cells / and j, and the notation j c n.n.(/) means 
that the inner summation is performed over the nearest 
neighbors of the /th cell only. We note in passing that 
this Heisenberg-like expression is valid only when the 
angles between neighboring moments are not too large; 



as shown in Ref. 17311 . neglecting this condition can lead 
to completely unphysical results. 

For a disordered system of finite elements having dif- 
ferent volumes, different distances between the element 
centers, and different numbers of nearest neighbors for 
each element, the expression Eq. (0 obviously can not 
be used. The most straightforward way to compute the 
exchange interaction in such a system would be to em- 
ploy a rigorous numerical approximation suitable for 
the evaluation of the integral Eq. (0, where the inte- 
grand values (magnetization vectors m) are given at ar- 
bitrarily placed spatial points. Derivation of such inte- 
gration formulas amounts to the approximation of the 
magnetization projections m v v z (r) using some kind of 
a polynomial interpolation of these functions between 
the points where their values are defined (in our case, 
polyhedron centers). The integrand in Eq. (0 includes 
first spatial derivatives of the magnetization field, so that 
corresponding polynomials for the energy evaluation 
should deliver a continuous (and even better smooth) 
first derivative of m(r). For the effective field evalua- 
tion, where the continuous expression involves second- 
order derivatives, the polynomial interpolation is even 
more demanding. In addition, we should keep in mind 
that the condition |m(r)| = 1 must be fulfilled every- 
where, so that the interpolation of the magnetization an- 
gles, rather than that of Cartesian magnetization compo- 
nents should be used. All these features would result in 
a highly complicated algorithm for the exchange energy 
and field evaluation, which might be, in addition, sub- 
ject to serious stability problems to usage of the angle 
interpolation. 

For these reasons, we have decided to develop an al- 
gorithm for the exchange-energy evaluation based on 
the summation of the nearest neighbors contributions 
[similar to the expression (Eq. |7]l] and at the same time 
take into account the differences between a regular grid 
and a disordered system mentioned above. To achieve 
this goal, we modify the expression Eq. (0 in the fol- 
lowing way. First, from the derivation of the expression 
Eq. © presented in 17311 . it is clear that for the regular 
mesh consisting of cubic cells, the volume AV in the nu- 
merator of this expression is actually not the volume of 
the cell, but the volume enclosed between the centers of 
the cell / and the neighboring cell j (for a cubic lattice or 
a lattice consisting of rectangular prisms, this volume is 
obviously equal to the cell volume, because it includes 
two halves of identical cells). Therefore, for a disor- 
dered system of arbitrary finite elements, AV should be 
replaced by Vy = (V, + V})/2 , where V, and Vj are the 
volumes of the /th and the jth finite elements. 
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The second adjustment of Eq. (J?) to a finite element 
system is the replacement of the distance a between the 
cell centers in a regular lattice by the distance Ar,- y be- 
tween the centers of cells i and j. 

The third and most complicated correction is due to 
the different number of nearest neighbors in a regular 
lattice and in a disordered system of finite elements. In 
a regular Cartesian lattice, each cell has exactly N m = 6 
nearest neighbors, and the angles between the lines con- 
necting the cell center with the centers of its neighbors 
in x, y and z directions are always 90°. For this rea- 
son, the overlapping of volumes enclosed between the 
centers of neighboring cells in, e.g., x and y directions 
is the same for all cells, what was taken into account 
by derivation of Eq. (O (it is important to note, that 
this overlapping should not be confused with the over- 
lapping of spheres mentioned in the discussion of algo- 
rithms used to contruct our disordered mesh in Sec. 12. U . 
In contrast to this nice feature of a regular Cartesian lat- 
tice, in a disordered system of finite elements the num- 
ber of nearest neighbors for different finite elements 
may be different, and the overlapping of volumes en- 
closed between the centers of a given cell and its differ- 
ent neighbors may also vary. For example, for the ele- 
ment with more than six nearest neighbors, the volumes 
enclosed between its center and centers of its neighbors 
would overlap more than for a cubic lattice. For such 
an element, the exchange-stiffness energy evaluated us- 
ing the sum Eq. would be overestimated due to this 
excessive overlapping, even when the two corrections 
explained above would be taken into account. 

The simplest method to solve this problem is the in- 
troduction of the correction factor 6/n av , where n av is 
the average number of nearest neighbors for the partic- 
ular random realization of our disordered finite-element 
system. This correction would compensate on aver- 
age the effect of the incorrect count of overlapping re- 
gions explained above. The accuracy of this simple 
correction method can be hardly estimated in advance, 
but both simple tests performed in [74] and additional 
much more complicated tests discussed in Sec. I2.3l be- 
low show that the accuracy provided by this correction 
method is surprisingly good. 

Summarizing, for magnetic moments belonging to 
the same phase we propose the following expression for 
the exchange-stiffness energy: 



1 N 

- z s 



2 Ar 2 . 

i=\ jcn.n.(i) i J 



2Aij Vjj , \ 
— ( m ' m /)' 



(8) 



where V,j = (V, + Vj)/2, Ary is the distance between 
the centers of the ith and the jth finite elements with 



volumes Vt and Vj, and Ay is the exchange constant. 

The last point to be discussed is the choice of nearest 
neighbors, which should be used in the inner summation 
in Eq. (0. The choice whether two elements should 
be considered as nearest neighbors is not unambiguous 
in disordered systems. We have adopted the following 
convention: two magnetic moments are considered as 
nearest neighbors, if the distance between the centers 
of corresponding polyhedrons is not larger than af max = 
1.4 (ji+rj). The cut-off factor / cut = 1 .4 is chosen so that 
for the overwhelming majority of finite elements those 
two of them which have a common face are treated as 
nearest neighbors. 

To evaluate the exchange-interaction energy between 
two finite elements (polyhedrons) belonging to different 
phases (hard and soft), we use the formula: 



J exch 



1 N 

!=1 jcn 



2A U V sp /2 
t & (Ar, 7 -fl hp ) : 



(m,m,). (9) 



Here, V sp is the volume of a soft-phase element and R^p 
is the radius of the sphere corresponding to the hard 
phase polyhedron. This modified expression Eq. ((9]) ac- 
counts for the fact that in this case the magnetization 
rotation occurs almost entirely within the polyhedron 
corresponding to the soft phase. 

2.2.3. Magnetodipolar interaction 

The energy of the long-range magnetodipolar interac- 
tion between magnetic moments and the corresponding 
contribution to the total effective field are computed us- 
ing the point-dipole approximation as 



£di P - -r YiftYi 



3e,j(eijltj)-lij 



i=i 



(10) 



i.e., magnetic moments of finite elements are treated as 
point dipoles located at the polyhedron centers. This ap- 
proximation is equivalent to the approximation of spher- 
ical dipoles, i.e., it would be exact for spherical finite el- 
ements. Hence, for our discretized system, this approx- 
imation introduces some computational errors, because 
our finite elements are polyhedrons. However, these er- 
rors are small, because the shape of these polyhedrons 
is close to spherical (see Fig. [3]), due to the special al- 
gorithm employed for the generation of our mesh, as 
explained in Sec. 12.11 

The summation in Eq. (TlOb is performed by the so 
called particle-mesh Ewald method. Didactically very 
instructive and detailed introduction into Ewald meth- 
ods can be found in Ref. [75]. The specific implemen- 
tation of the lattice-based Ewald method for the magne- 
todipolar interaction for regular and disordered systems 
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of magnetic particles is described in our papers 

Here, we briefly repeat the basic issues of this algorithm 

to make our paper self-containing. 

First we remind that the Ewald method IT78ll was ini- 
tially invented for evaluating conditionally converging 
lattice sums for the Coulomb interaction in ionic crys- 
tals. At present, it is a standard method to calculate 
any long-range interaction — including Coulomb sums, 
gravitation energy, dipole interaction, elastic forces in 
dislocation networks etc. — in systems with periodic 
boundary conditions (PBC). In such systems direct sum- 
mation over all field sources is impossible simply due 
to their infinite number. Hence, we must use a Fourier 
expansion over the reciprocal lattice vectors k which 
correspond to the infinitely repeated simulation volume. 
For the point sources of the long-range field, the cor- 
responding Fourier components decay relatively slowly 
with increasing magnitude of the wave vector k in recip- 
rocal space. In numerical simulations we always have 
to our disposal only a finite number of such wave vec- 
tors, so that the Fourier spectrum of our long-range in- 
teraction should be cut off at the maximal finite value 
kmax- As mentioned above, the Fourier harmonics de- 
cay slowly, so that at k m . dx they are by no means small. 
For this reason, the spectrum cut-off due to the elimina- 
tion of all Fourier components with k > £ max is sharp, 
thus, leading to large artificial oscillations of the interac- 
tion potential after its inverse transformation to the real 
space. 

As with nearly all Ewald methods, the version de- 
scribed below for dipolar systems solves the problem 
by adding and subtracting a Gaussian dipole at each lo- 
cation of a point dipole in the initial system. Using 
the definition of the gradient of the 6 function, it is easy 
to show that this operation corresponds to the addition 
and subtraction of a charge distribution (with a width <x) 



P?(r) 



exp 



(r - rjT 

2cr 2 



(11) 



H* p + H* p is 



Then the magnetodipolar field H dlp 
evaluated as the sum of two contributions from subsys- 
tems A and B. The first subsystem consists of Gaussian 
dipoles Eq. (fTTT i and the second one is composed of the 
original point dipoles minus these Gaussian dipoles, 



PB(r) = -Y j [ f i i V8(r-r i )-pf(r)] ! 



(12) 



where the first terms in the sum on the right represent 
the charge density of a point dipole located at r,-. The 



field created by a composite object in square brackets of 
Eq. (HI is [76] 



#B,/( r - r '-) 



3(a-q i )(n i Ar i ) p.? 
Arf Ar 3 

2 (a - ad (fi; Ar,) 



Ar 5 



exp 



/o(Ar,) + 

Ar 2 
'la 2 



(13) 



where a — x, y, z. 

It is important to note that the function foif) decays 
with distance as exp(-r 2 ), 



Mr) = 1 - erf 



■V2 



V2r 
+ exp 



r 

'la 2 



(14) 



The goal of this decomposition of the original system 
of point dipoles is the following. The field Eq. ( IT3b from 
the second subsystem B is a short-range one, because 
each point dipole is screened by a Gaussian dipole with 
the same total moment, but with the opposite sign. The 
computation of such a short-range contribution takes 
~ N operations for a system of N particles. The first 
subsystem A consists of dipoles having a smooth Gaus- 
sian charge distribution, so that its Fourier components 
decay rapidly with increasing k. This fast decay allows 
a painless cut-off of the Fourier spectrum at large wave 
vectors, so that the contribution from the first subsystem 
Eq. (fTTT i can be safely calculated using Fourier expan- 
sion. More detailed explanations concerning this proce- 
dure can be found in Ref. [76]. 

Already the above most straightforward implementa- 
tion of the Ewald method allows for a reliable evalua- 
tion of the dipolar field in systems with PBC. However, 
for disordered systems, this method has the same pro- 
hibitively high operation count ~ N 2 , as a direct summa- 
tion for systems with OBC. The reason is that particle 
positions in disordered systems do not form a regular 
lattice, so that the Fourier transformation for the cal- 
culation of the long-range contribution H* p can not be 
done via the fast Fourier transformation technique: ex- 
ponential factors exp(ikr,) should be computed for all 
wave vectors k and all particle positions r,- separately, 
leading to the operation count given above. 

In order to decrease the computational costs, several 
lattice versions of the Ewald method have been devel- 
oped (see the overview H)- The general idea behind 
all these methods is to employ some mapping of the ini- 
tial disordered system onto a regular lattice, in order 
to enable the application of the FFT Using this gen- 
eral paradigm, we have implemented the following al- 
gorithm: 
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(z) First, we map our disordered system of point mag- 
netic dipoles fii = //(r,) onto a system of dipoles lo- 
cated at lattice points r p (p is the 3D index) using some 
weighting function W3</(r), 

N 

#0" P ) = ^J"(ri)w M (|r; - r p \) = 
1=1 

2^Hi w(\Xi - x p \) wQyt - y p \) w(\zj - z p \). (15) 

1=1 

We emphasize that the whole method makes only sense 
if the mapping function w is strongly localized, so that 
the sum over all N dipoles in Eq. $15[ is actually re- 
stricted to a few nearest neighbors M„b of the lattice 
node p. 

(z'z) Next, we add and subtract to each point lattice 
dipole two Gaussian dipoles Eq. d!2t . as in the straight- 
forward Ewald method described above. 

(z'zz) Further, we compute the dipolar field of this lat- 
tice system as described above, i.e., as the sum of the 
long-range contribution from smooth Gaussian dipoles 
positioned on the lattice and the short-range contribu- 
tion Eq. ( fT3l from the composite objects "point dipole 
- Gaussian dipole", also placed on the lattice. 

(z'v) Finally, the field obtained in this way on the lat- 
tice points r p is mapped back onto the initial dipole lo- 
cations r, using the same functions w as in Eq. (TT3T >. 

As mentioned above, the major advantage of this lat- 
tice Ewald version is the possibility to use FFT for com- 
puting the long-range part of the total magnetodipolar 
field. In addition, we can also accelerate the evalutaion 
of the short-range contribution. Namely, we note that (/) 
the contribution Eq. ( ITJb depends only on the difference 
Ar between the source and target coordinates and (if) 
both source and target points are located on the lattice. 
Hence, this short-range contribution also is a discrete 
convolution and as such can be also computed by the 
FFT technique. Using this nice feature, we can increase 
the number of nearest neighbor shells used by the eval- 
uation of the short-range interaction part without addi- 
tional time cost, making the corresponding truncation 
error arbitrarily small. Keeping in mind that the evalu- 
ation of the long-range field part via the FFT technique 
for the lattice system is exact, we conclude that the only 
source of computational errors in our algorithm is the 
mapping of the initial disordered system onto a lattice, 
which can be easily controlled and reduced by choos- 
ing the suitable mapping scheme IT79ll . We have found 
that already the conventional first-order mapping used 
together with a lattice having a cell size equal to R sp /2 
(here /? sp is the sphere radius used to generate the mesh 



for the soft-phase discretization) ensures by the evalua- 
tion of H dlp a relative error smaller than 0.01, which is 
good enough for our purposes. 

2.3. New micromagnetic algorithm: minimization pro- 
cedure and numerical tests 
For the minimization of the total magnetic energy, ob- 
tained as the sum of all contributions described above, 
we use the simplified version of a gradient method 
employing the dissipation part of the Landau-Lifshitz 
equation of motion for magnetic moments [65, 80]. This 
means that we update the magnetization configuration at 
each step as 

m new = m old _ At ^old x ^old x jjeffjj > (lfi) 

where m,- denotes the unit magnetization vector m, = 
M,- 1 Ms and h eff is the reduced effective field, evaluated 
in a standard way as the negative energy derivative over 
the magnetic moment projections ll65ll . 

Since we are looking for the energy minimum, the 
time step in Eq. ( TToT l is chosen and adapted using the 
monitoring of this energy. If the total energy decreases 
after the iteration step performed according to Eq. ( fTSI l. 
we accept this step. If the energy increases, we restore 
the previous magnetization state, halve the time step 
(Af — > Af/2) and repeat the iteration. During the mini- 
mization procedure we may also increase the time step 
to avoid an unnecessary slow minimization: the time 
step is doubled, if the last few steps (typically 5-10 
steps) were successful. For the termination of the min- 
imization procedure, we use the local torque criterion: 
we stop the iteration process, if the maximal torque act- 
ing on magnetic moments is smaller than some pre- 
scribed value, i.e., max,-{|m,- x h eff |) < s. As is well 
known, this condition is more appropriate than the alter- 
native criterion of a sufficiently small energy difference 
between the two subsequent steps. In all tested cases 
the value s = 10~ 3 was small enough to ensure the min- 
imization convergence. 

The new methodology explained in detail above was 
tested on two simple examples in Ref . IT74fl . We remind 
that we have first reproduced — using our disordered 
mesh — with a high accuracy the analytically known 
magnetization profile of a standard 3D Bloch wall. Sec- 
ond, for a trial 3D magnetic configuration defined via 
simple trigonometric functions of coordinates (we have 
used these functions to ensure a sufficiently slow spa- 
tial variation of the system magnetization direction), we 
have obtained a very good agreement between the to- 
tal energy and partial energy contributions found by our 
new method and the FDM micromagnetic package Mi- 
croMagus lHHl 
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Figure 2: (a) Vortex (2D cross section) and (b) flower (3D arrow plot) magnetization states obtained by our new methodology explained in Sec. 12. 31 
above. Data taken from Ref. 1 821]. 



Here, we would like to present two additional much 
more complicated tests, where we compute the equilib- 
rium magnetization configurations of a cubic magnetic 
particle obtained using our new method and compare 
these configurations with the results obtained for the 
same system by the MicroMagus package. 

The particle size was chosen to be 40 x 40 x 40 nm 3 , 
and the magnetic materials parameters were set to Ms = 
800kA/m,A = l.Ox 10" 11 J/m, and K = 5.0x 10 4 J/m 3 
(uniaxial anisotropy). For the simulations using our new 
method, the particle was discretized into N = 9000 
polyhedrons with a typical size of d = 2 nm. For the 
MicroMagus simulations, a cubic cell with a side length 
of 2.5 nm was used. 

For the test problems, we have chosen two well 
known magnetization states typical for ferromagnetic 
particles of this size l83ll : the vortex state and the so- 
called flower state. To obtain the vortex state, we have 
started the minimization procedure from the state that 
is topologically equivalent to the vortex, the so-called 
closed Landau domain configuration. The flower state 
could be obtained by starting the energy minimization 
simply from the homogeneous configuration with the 
magnetization directed along a cube side. 

Table Q] lists the total energies, partial energy contri- 
butions and the reduced magnetization values for the 
equilibrium magnetization states shown in Fig. [2] ob- 
tained by the new method and by the standard FDM 
simulations (MicroMagus package). Almost all energy 
contributions obtained by the two methods agree very 
well. The only significant relative difference can be 
found for the anisotropy energy of the flower state; how- 
ever, this significant relative difference {AE/E) arises 



Table 1 : Comparison of energies and reduced magnetizations for the 
vortex and flower magnetization states computed by the new method 
and by the standard finite difference simulations (MicroMagus soft- 
ware). Data taken from Ref. 1 82J|. 



Vortex energies (xlO 


l!i J) 


New method 


MicroMagus 


Etot 




8.225 


8.270 


^an 




1.361 


1.385 


^exch 




4.409 


4.562 






2.455 


2.324 


M/M s 




0.400 


0.406 




Flower energies (xlO 


1S J) 


New method 


MicroMagus 


Etot 




7.813 


7.843 


E;m 




0.137 


0.127 


Eexch 




0.434 


0.441 


Edip 




7.242 


7.275 


M/M s 




0.972 


0.974 



simply due to a very low value of this energy. All in 
one, the agreement between the new and the established 
methodologies for all cases where the standard methods 
are applicable is fully satisfactory. 

3. Magnetic SANS cross section of unpolarized neu- 
trons 

In our micromagnetic simulations of elastic magnetic 
SANS we have focussed on the two most commonly 
employed scattering geometries where the wavevector 
ko of the incident neutron beam is either perpendicu- 
lar [case (i)] or parallel [case (ii)] to the external mag- 
netic field H, which is applied along the z-direction of 
a Cartesian coordinate system. (e. v , e y and e z repre- 
sent the unit vectors along the Cartesian axes.) Further- 
more, since the focus of the present study is on mag- 
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Figure 3: Sketch of the two scattering geometries and of the micro- 
scopic structure of the nanocomposite sample (simulation volume: 
V = 250 X 600 X 600 nm 3 ). Blue polyhedrons — Fe particles; yellow- 
orange-red polyhedrons — matrix phase. In the micromagnetic simu- 
lations the sizes of the Fe particles are distributed according to a Gaus- 
sian function centered at about 10 nm and the particle volume fraction 
equals 40 %, as in the experimental study Ref. (1711 . 

netic spin-misalignment scattering, we have ignored 
the nuclear SANS contribution. Note, however, that 
for polycrystalline texture-free magnetic nanocompos- 
ites the nuclear SANS signal is virtually independent of 
the applied magnetic field and isotropic, and its mag- 
nitude is generally small compared to the here relevant 
spin-misalignment scattering [31]. Furthermore, the re- 
striction to unpolarized neutrons entails the neglect of 
scattering contributions from helical spin arrangments, 
which are of relevance, e.g., in FeCoSi and MnSi single 
crystals lEjEj]. 

A sketch of the above two scattering geometries along 
with a schematic drawing of the microscopic structure 
of the nanocomposite sample can be seen in Fig. [3] 

3.1. Case (i): k ± H || e. 

For ko || e x , the elastic magnetic SANS cross section 
dI,M/dCl at momentum- transfer vector q reads 113 111 

^(q) = ^4(|M,| 2 + |M v | 2 cos 2 e 
+\M Z \ 2 sin 2 9 - (M y M* z + M*M Z ) sin 9 cos 6>), (17) 



where V = 250 x 600 x 600 nm is the scattering vol- 
ume, bu = 2.9 x 10 8 A^'irT 1 , c* is a quantity complex- 
conjugated to c, and M( WlZ )(q) are the Fourier trans- 
forms of the magnetization components M( A V Z )(r). Note 
that in the small-angle limit and for this particular ge- 
ometry the scattering vector q can be expressed as q = 
q (0, sin 9, cos 9), where 9 denotes the angle between q 
and H (Fig. H. 

3.2. Case (ii): k || H || e z 
For this geometry, one finds 

^(q) = ^ 2 (|M,| 2 sin 2 + |M r | 2 cos 2 

+\M Z \ 2 - (M x My + M* x My) sin 9 cos 9), (18) 

where q = q (cos 9, sin 9, 0) and 9 is measured relative 
toe, (Fig. [3). 

For the micromagnetic simulations of magnetic 
SANS from two-phase nanocomposites, we used (un- 
less otherwise stated) the following materials param- 
eters for hard ("h") and soft ("s") phases, which are 
characteristic for the Fe-based nanocrystalline alloy 
NANOPERM H: magnetizations M h = 1750 kA/m 
and M s = 550 kA/m, anisotropy constants = 4.6 X 
10 4 J/m 3 and K s = 1.0 X 10 2 J/m 3 . As a value for the 
exchange-stiffness constant we used A = 0.2xl0~ n J/m 
for interactions both within the soft phase and between 
the hard and soft phases. 

4. Results and discussion 

The applied-field dependence of the total magnetic 
SANS cross sections dI,MldQ. [computed, respectively, 
by means of Eqs. ( fTTT i and (fT8l l is displayed in Fig.|4]for 
both scattering geometries, i.e., for the situations when 
the wavevector ko of the incoming neutron beam is per- 
pendicular [case (i)] or parallel [case (ii)] to the applied 
magnetic field H, which for both cases is assumed to be 
parallel to e z . The corresponding radially-averaged data 
can be seen in Fig. [5] 

While dT.M/dQ. for ko || H is isotropic (i.e., 9 inde- 
pendent) over the whole field and momentum-transfer 
range, it is highly anisotropic for ko _L H (Fig. At 
a saturating applied magnetic field of fioH - 1 .5 T, 
where the normalized "sample" magnetization is (for 
both geometries) larger than 99.9 %, the anisotropy of 
dl,M/dD. [case (i)] is clearly of the sin 2 0-type, i.e., elon- 
gated normal to H; this is because magnetic scattering 
due to transversal spin misalignment is small close to 
saturation and the dominating scattering contrast arises 



1500 mT 290 mT lOOmT 30 mT 




Figure 4: Applied-field dependence of the total magnetic SANS cross section dLu/dSl for ko iH (Eq. )17t , upper row) and for ko || H (Eq. 1 1 81 . 
lower row). The external magnetic field H || e. is applied horizontally in the plane of the detector for ko ± H (q x = 0, upper row) and normal to 
the detector plane for ko || H (q z = 0, lower row). Materials parameters of NANOPERM were used (see text). Pixels in the corners of the images 
have q = 1.2 mrr 1 . Logarithmic color scale is used. 




Figure 5: Solid lines: Radially-averaged dZu/dil as a function of scattering vector q for ko X H (left image) and for ko || H (right image) (data 
have been smoothed). Field values (in mT) from top to bottom, respectively: 30, 100, 290, 1500. Dashed lines in both images: dHu/dil oc q~ 4 . 
Solid circles in both images represent part of the (squared) form factor of a sphere with a radius of R = 5.7 nm. 



from nanoscale jumps of the longitudinal magnetiza- 
tion at phase boundaries. On decreasing the field, the 
transversal magnetization components increase in mag- 
nitude as long-range spin misalignment develops at the 
smallest q. The SANS pattern in case (i) essentially re- 
mains of the sin 2 f9-type at lower fields, although a more 
complicated anisotropy builds up at small q. As can be 
seen in Fig. [5J dl.M/dQ. at small q increases by more 
than one order of magnitude as the field is decreased 
from 1.5 T to 30 mT. Asymptotically, at large q, the 
power-law dependence of dl.M/dQ. can be described by 



dl,M/d£l oc q~ 4 . In agreement with the nature of the un- 
derlying microstructure (~ lOnm-sized single-domain 
Fe particles in a nearly saturated matrix), one can de- 
scribe the oscillations of dY-M/dQ. at large q and H by 
the form factor of a sphere with a radius of R = 5.7 nm 
(solid circles in Fig. [5]). We remind that the shape of the 
particles in our micromagnetic algorithm is not strictly 
spherical. 

Figure [6] shows the projections of the magnetization 
Fourier coefficients \M X \ 2 , |M,,| 2 , \M Z \ 2 , and of the cross 
terms CT = -(M,M* + M*,M Z ) and CT = -(M X M* + 
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Figure 6: Results of the micromagnetic simulations for the Fourier coefficients of the magnetization. The images represent projections of the 
respective functions into the plane of the 2D detector, i.e., q x = for ko 1 H (left image) and q z = for ko || H (right image). From left column to 
right column, respectively: |M V | 2 , |M } ,| 2 , \M Z \ 2 , and CT = -(M y Mt + M*M.) (left image) and CT = -(M,-M* + M*M y ) (right image). In the first 
three columns from left, red color corresponds, respectively, to "high intensity" and blue color to "low intensity"; in the fourth column, blue color 
corresponds to negative and orange-yellow color to positive values of the CT. Pixels in the corners of the images have q = 1.2 nm -1 . Logarithmic 
color scale is used. 

M* x My) into the plane of the 2D detector at the same 
external-field values as in Fig. [4] It can be seen in 
Fig.|6]that in case (i) both \M X \ 2 and \M Z \ 2 are isotropic 
over the displayed (q, H) range, while the Fourier coef- 
ficient \M V \ 2 reveals a pronounced angular anisotropy, 
with maxima that lie roughly along the diagonals of 
the detector (the so-called "clover-leaf" anisotropy, see 
Fig. Qj] below). In case (ii), \M X \ 2 and |M,,| 2 are both 
strongly anisotropic (with characteristic maxima in the 
plane perpendicular to H), while \M Z \ 2 is isotropic. 
When [for case (ii)] all Fourier coefficients are multi- 
plied by the corresponding trigonometric functions and 
summed up [compare Eq. (1181) 1. the resulting dLuldSl 
becomes isotropic (Fig.|4j lower row). 

The cross terms for both scattering geometries vary 
in sign between quadrants on the detector. The respec- 
tive CT is positive in the upper right quadrant of the 
detector (0° < 9 < 90°), negative in the upper left quad- 
rant (90° < 9 < 180°), and so on. When both CT's 
are multiplied by sin 9 cos 9, the corresponding contri- 
bution to dl-M/dQ. becomes positive-definite for all an- 
gles 9 (compare Fig. |7j. This obervation suggests that — 
contrary to the common assumption that the CT av- 
erages to zero for statistically isotropic polycrystalline 
microstructures — the CT appears to be of special rele- 



vance in nanocomposite magnets. 




Figure 7: (left image) CT = -(M y M* + M*,M Z ) for k ± H (p Q H = 
0.29 T). Blue color corresponds to negative and orange-yellow color 
to positive values of the CT. (right image) CT sin 9 cos 0. Red color 
corresponds to "high intensity" and blue color to "low intensity". All 
other settings are as in Fig. [6] 

In Fig. [8] we show for both scattering geometries 
the radially-averaged total dEju/dSl along with the 
radially-averaged individual scattering contributions to 
dLni/dQ, i.e., the radial average of terms ^■bj 1 \M x \ 2 , 
^-b 2 H \M y \ 2 cos 2 9, ^-b 2 H CT sinfl sin0, and so on [com- 
pare Eqs. (Ell and (18])]. At saturation (p. H = 1.5 T), 
both transversal scattering contributions, i.e., terms cc 
\M x (q)\ 2 and oc \M y (q)\ 2 , are for both cases (i) and (ii) 
small relative to the other terms and the main contribu- 
tion to the total dl.M/dQ. originates from longitudinal 
magnetization fluctuations, i.e., from terms oc \M z (q)\ 2 . 
For ko || H, both transversal terms are so small that they 



12 



k Q _LH; 1500 mT 




q (nm ) 



k _L H; 290 mT 




q (nm ) 
k Q 1 H; 30 mT 






io 4 






1 

s 

u 


io 3 


g 


io 2 




1 



k Q II H; 1500 mT 



g 



10' 



0.1 0.2 



0.7 



q (nm ) 
k Q II H; 290 mT 




q (nm ) 
k Q II H; 30 mT 



q (nm ) 




q (nm ) 



Figure 8: Radially-averaged total magnetic SANS cross sections dZ^/dil and radially-averaged individual contributions to dZu/dQ. as a function 
of scattering vector q and applied magnetic field H for ko 1 H and ko || H (see insets). Black lines: total dEu/dil; Blue lines: |M,| 2 ; Magenta 
lines: CT; Green lines: |M,,| 2 ; Red lines: \M X \ 2 . In the above notation, the prefactor ^-b 2 H and the numerical factors that result from the averaging 
procedure have been omitted for clarity. Note that different trigonometric functions may be involved in the averaging procedure, compare, e.g., 
^b 2 H \M x \ 2 for k ± H [Eq. (TTJ] and ^-b 2 H \M x \ 2 sin 2 9 for ko || H [Eq. flStl. 



are not visible within the displayed "intensity" range 
and dl.M/dQ. practically equals the \M z (q)\ 2 scattering 
(both curves superimpose). Note that the CT for case (i) 
is the product of a transversal and the longitudinal mag- 
netization Fourier coefficient, whereas for case (ii) the 
CT contains the two transversal components. This ex- 
plains why the CT for case (ii) is much smaller than the 
CT for case (i) at fields close to saturation. On decreas- 
ing the field, the transversal Fourier coefficients and the 
CT's become progressively more important, in particu- 
lar at small q. 



It is also important to note that the present simula- 
tions were carried out by assuming a quite large jump 
in the magnetization magnitude AM at the interphase 
between particles and matrix, /uoAM = 1 .5 T. Conse- 
quently, the ensuing \M z (q)\ 2 scattering in both geome- 
tries and the CT scattering in case (i) are relatively large. 
For AM = 0, the CT at saturation for case (i) becomes 
negligible, since M z (q) oc d(q = 0). 

Figures [6] and [8] embody the power of our approach: 
By employing numerical micromagnetics for the com- 
putation of magnetic SANS cross sections, it becomes 
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possible to study the individual magnetization Fourier 
coefficients and their contribution to dl.M/dCl. This 
sheds light on the ongoing discussion regarding the 
explicit q-dependence of dLuid£i 15611 . In particu- 
lar, the approach of combining micromagnetics and 
SANS complements neutron experiments, which gener- 
ally provide only a weighted sum of Fourier coefficients 
[compare Eqs. ( fT77T > and ([T811. a fact that often hampers 
the straightforward interpretation of recorded SANS 
data. While it is in principle possible to determine some 
Fourier coefficients, e.g., through the application of a 
saturating magnetic field or by exploiting the neutron- 
polarization degree of freedom via so-called SANSPOL 
or POLARIS methods (e.g., Refs. |EE,Ellls3l), it is dif- 
ficult to unambiguously determine a particular scatter- 
ing contribution without "contamination" by unwanted 
Fourier components. For instance, when the applied 
field is not large enough to completely saturate the sam- 
ple, then the scattering of unpolarized neutrons along 
the field direction does not represent the pure nuclear 
SANS, but contains also the magnetic SANS due to the 
misaligned spins 03911 . 



290 mT 



100 mT 



30 mT 



The finding [for case (i)] that \M X \ 2 and \M : \ 2 are 
isotropic and that |M V | 2 = \M v \ 2 (q, ff) provides a straight- 
forward explanation for the experimental observation 
of the clover-leaf anisotropy in the SANS data of 
the nanocrystalline two-phase alloy NANOPERM lfl7ll . 
Our simulation results for the difference cross section 
oc (\M X \ 2 + \M y \ 2 cos 2 6 + CT sin 6 cos 0) (see Figs.|9]and 
[Tol l, where the scattering at saturation (fXoH = 1 .5 T) has 
been subtracted, agree qualitatively well with the ex- 
perimental data 11741 18211 . Clover-leaf-type anisotropies 
in dI.M/dCl have also been reported for a number of 
other materials, including precipitates in steels 1391. 
nanocrystalline Gd |32l |35J, and nanoporous Fe 113 311 . 



The maxima in the difference cross section [for 
case (i)] depend on q and H, and may appear at angles 6 
significantly smaller than 45°. This becomes evident in 
Fig. Q~TJ where we show (for ko ± H) polar plots of the 
simulated difference cross section at selected q and H. 

The results of our previous work 11821 18511 strongly 
suggest that the magnetodipolar interaction plays a de- 
cisive role for the understanding of magnetic SANS of 
nanocomposites. In fact, it is this interaction which is 
responsible for the anisotropy, i.e., for the 0-dependence 
of the magnetization Fourier coefficients and, hence, 
of dI.M/dQ.. The impact of the dipolar interaction on 
dI.M/dQ. can be conveniently studied, since our micro- 
magnetic algorithm allows one to "switch on" and "off" 
this energy term. Figure [12] shows results of micro- 
magnetic simulations for \M y \ 2 and for both CT's ob- 




Figure 9: Comparison between simulation (upper row) and experi- 
mental data (lower row) for the difference cross section oc (|A/ r | _ + 



\MyY- 



- 9 + CT sin 6 cos 8) at different external fields as indicated 



(ko -L H). Pixels in the corners of the images have q = 0.64 nm~'. 
Logarithmic color scale is used. Since the experimental data was not 
obtained in absolute units, we have multiplied it with a scaling factor 
for comparison with the simulated data. H is horizontal in the plane. 
Experimental data were taken from Ref. fl7ll . 
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Figure 10: (•) Radially-averaged experimental difference cross sec- 
tions as a function of momentum transfer q and H (ko ± H). Field 
values (in mT) from top to bottom: 30, 100, 290. Solid lines: Results 
of the micromagnetic simulations (data have been smoothed). Verti- 
cal dashed lines indicate the region where the clover-leaf anisotropy is 
observed. Experimental data were taken from Ref. [17] and multiplied 
by a scaling factor (compare Fig.[5J- 



tained with and without the dipolar interaction. When 
the dipolar interaction is ignored in the micromagnetic 
computations, all Fourier coefficients are isotropic at all 
q and H investigated (data for \M X \ 2 and \M Z \ 2 are not 
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Figure 11: Polar plots of the simulated difference cross section oc 
(|M V | 2 + |M V | 2 cos 2 + CT sin 6 cos 9) at different combinations of mo- 
mentum transfer q and applied magnetic field H (see insets) (ko ± H). 
Data have been smoothed. Dotted lines (oc sin 2 cos 2 9) serve as 
guides to the eyes. 

shown). This observation shows that for any realistic 
description of experimental magnetic SANS data this 
interaction has to be taken into account. 

Generally, the sources of the magnetodipolar field 
are nonzero divergences of the magnetization (V • M + 
0). For magnetic nanocomposites, the most prominent 
"magnetic volume charges" are related to the nanoscale 
variations in the magnetic materials parameters at the 
phase boundary between particles and matrix, e.g., vari- 
ations in the magnetization, anisotropy or exchange in- 
teraction. Such jumps in the magnetic materials pa- 
rameters may give rise to an inhomogeneous spin struc- 
ture which decorates each nanoparticle. Figure [T3] dis- 
plays the real-space magnetization distribution around 
two nanoparticles. Note that the symmetry of the spin 
structure replicates the symmetry of the CT (compare to 
Fig HJ. In the presence of an applied magnetic field the 
stray-field and associated magnetization configuration 
around each nanoparticle "look" similar (on the aver- 
age), thus giving rise to dipolar correlations which add 
up to a positive-definite CT contribution to dY^MldQ.. 
Note, however, that for polycrystalline microstructures 
clover-leaf-type anisotropies may become only visible 
in dl-nt/dCl for ko _L H. 

As mentioned above, not only variations in the mag- 
netization magnitude, but also variations in the direc- 
tion and/or magnitude of the magnetic anisotropy K 
(random anisotropy) and variations in the magnitude of 
the exchange coupling may give rise to dipolar corre- 
lations. The micromagnetic simulation package allows 



M I 2 CT 




Figure 12: Influence of the dipolar interaction on the Fourier co- 
efficients of the magnetization. |M V | 2 and both cross terms CT = 
-(M y M* z + M*M Z ) (k X H) and CT = -(M X M* + M*M y ) (k || H) 
were computed from a real-space magnetic microstructure with a nor- 
malized magnetization of 99.0%. Applied fields of 290 mT (with 
dipolar interaction) and 7 mT (without dipolar interaction En p = 0) 
were required in order to achieve this magnetization value. The cor- 
responding results for the Fourier coefficient \M X \" (for ko || H) are 
analogous to the depicted results for |M,,| 2 . Pixels in the corners of the 
images have q = 0.9 nm -1 . Logarithmic color scale is used. 

us to vary the magnetic materials parameters of both 
phases of the nanocomposite. Hence, it becomes pos- 
sible to study the impact of such situations on the mag- 
netic SANS. 

In order to investigate variations in K (which are, by 
construction, naturally included into our micromagnetic 
algorithm), we have computed the spin distribution for 
the situation that M h = M s = M (i.e., AM = 0), but for 
different values of M. Figure [14] reveals that a clover- 
leaf-type pattern in |M,.| 2 develops with increasing mag- 
netization value M, i.e., with increasing strength of the 
magnetodipolar interaction. As jumps in the magneti- 
zation at phase boundaries are excluded here as possi- 
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Figure 13: Results of a micromagnetic simulation for the 2D spin distribution around two selected nanoparticles (blue circles), which are assumed 
to be in a single-domain state. The external magnetic field H is applied horizontally in the plane (poH = 0.3 T). Left image: Magnetization 
distribution in both phases; note that /ioAM = /io(Mh _ M s ) = 1.5 T. In order to highlight the spin misalignment in the soft phase, the right image 
displays the magnetization component M ± perpendicular to H (red arrows). Thickness of arrows is proportional to the magnitude of Mj_. Blue 
lines: Dipolar field distribution. 



ble sources for perturbations in the spin structure, it is 
straightforward to conclude that nanoscale fluctuations 
in K give rise to inhomogeneous magnetization states, 
which decorate each nanoparticle and which look sim- 
ilar to the structure shown in Fig. [13] This observation 
strongly suggests that the origin of the clover-leaf pat- 
tern in (TLbi/dO. of nanomagnets is not only related to 
variations in magnetization magnitude but also due to 
variations in the magnitude and direction of the mag- 
netic anisotropy field. 

250kA/m 500 kA/m 750 kA/m 1030 kA/m 




Figure 14: Fourier coefficient \M y \ (q) at fioH = 0.3 T and for = 
M s = M (i.e., AM = 0) (kn ± H). M increases from left to right 
(see insets). K h = 4.6 X 10 4 J/m 3 , K s = 1.0 X 10 2 J/m 3 and random 
variations in easy-axis directions from particle to particle are assumed. 
Data taken from Ref. (H. 



5. Summary and conclusions 

By means of a recently developed micromagnetic 
simulation methodology — especially suited for model- 
ing multi-phase materials — we have computed the mag- 
netic small-angle neutron scattering (SANS) cross sec- 
tion dliM/dCl of a two-phase nanocomposite magnet 
from the NANOPERM family of alloys. Besides taking 
into account the full nonlinearity of Brown's equations 



of micromagnetics, the approach allows one to study 
the dependency of the individual magnetization Fourier 
coefficients M( X ,y, Z ) on the applied magnetic field H 
and, most importantly, on the momentum-transfer vec- 
tor q. This ideally complements neutron experiments, in 
which a weighted sum of the M( X ,y,z) lS generally mea- 
sured. It is this particular circumstance, in conjunc- 
tion with the flexibility of our micromagnetic package 
in terms of microstructure variation (particle size and 
distribution, materials parameters, texture, etc.), which 
makes us believe that the approach of combining full- 
scale 3D micromagnetic simulations with experimental 
magnetic-field-dependent SANS data will provide fun- 
damental insights into the magnetic SANS of a wide 
range of magnetic materials. The micromagnetic simu- 
lations underline the importance of the magnetodipolar 
interaction for understanding magnetic SANS. In partic- 
ular, the so-called clover-leaf-shaped angular anisotropy 
in dZtf/dQ — which was previously believed to be ex- 
clusively related to nanoscale jumps in the magnetiza- 
tion magnitude at internal interphases — is of relevance 
for all bulk nanomagnets with spatially fluctuating mag- 
netic parameters. 
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